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Molecular constants have been computed 
for the ground states ^U of ^^OH and ^X 
of ^°^AgH^. The valence-bond method and 
advanced computational technique were 
used to perform all-electron ab initio cal- 
culation of molecular electronic struc- 
tures. The basic idea behind the model is 
to introduce the molecular wave functions 
in terms of Hartree-Fock many-electron 
atomic determinants. Full configuration 
interaction (CI) with nonorthogonal basis 
leads to the accurate calculation of 
molecular constants such as dissociation 
energy, equilibrium bond distance, vibra- 
tional and rotational constants with an 



agreement to the experimental data 
within a few percent. 



Key words: AgH^ molecule; diatomic 
molecule; dissociation energy; equilibrium 
bond distance; molecular constants; OH 
molecule; valence-bond method; vibrational 
and rotational constants. 



Accepted: February 20, 1998 



1. Introduction 



We present a new version of the valence-bond (VB) 
method, which is able to calculate the electronic 
structure of diatomic molecules. The first version of this 
ab initio method was developed by Heitler and London 
[1] and applied to the H2 molecule. The main idea of the 
method is to describe the molecular wave function as a 
linear combination of localized atomic wave functions. 
The new version of the method is oriented toward the 
calculation of many-electron molecules [2]. The self- 
consistent wave functions of the constituent atoms are 
constructed from determinants by using the Hartree- 
Fock (HF) method for non-relativistic atoms or the 
Dirac-Fock (DF) method for relativistic atoms. The con- 
struction of molecular wave functions from the many- 
electron atomic wave functions leads to the proper disso- 
ciation behavior. The total energy of the molecule as a 
function of inter-nuclear separation is presented as a 



power series in overlap integrals. This series converges 
quickly with atomic HF and DF functions. 

The molecular basis set contains both occupied and 
unoccupied atomic HF or DF orbitals. Since the excited 
orbitals have a large radius and are not effective in 
describing the correlation and polarization in the 
molecule, we create additional basis functions made 
of localized orbitals, i.e., the Sturm's wave functions 
[3]. This way of constructing the molecular wave 
functions avoids the need for large basis sets and ensures 
much shorter computation times as well as smaller 
errors in the calculation of the correlation expansion. 
The configuration interaction (CI) procedure, which 
takes into account the most important inner and 
outer shell correlation, simplifies and converges 
quickly as well, due to presence of HF or DF occupied 
orbitals. 
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The important advantage of this method is that the 
atomic core orbitals are not introduced by a pseudo- 
potential but calculated exactly. Although the inner-core 
orbitals do not affect significantly the molecular forma- 
tion, the exact representation of the core orbitals is 
important for properties of molecules that depend on the 
inner core wave functions such as, for example, the 
hyperfine interactions. These core effects are especially 
pronounced for molecules with heavy atoms. In princi- 
ple the Coulomb and exchange interactions can be 
calculated between all atomic orbitals, but this proce- 
dure can be limited to the most important orbitals to 
simplify the diagonalization procedure and save 
computer memory for the CI. The remaining orbitals 
will then form the core. A crucial step in the method 
is the calculation of many-center integrals using the 
sophisticated HF and DF functions. For the dimers the 
evaluation of the two-center integrals is simplified by 
using a reexpansion procedure [4]. In this procedure 
each of the two-center integrals is expanded in terms of 
one-center integrals located on the atomic centers. We 
have found that convergence of this expansion is acceler- 
ated by dividing the range of integration in two and 
considering only a reexpansion of the slowly varying 
part of the atomic wave function away from its nuclear 
center. 

The formalism of the VB method is described in 
more detail in Ref. [2]. Here we discuss the computa- 
tional details of a calculation of the molecular constants 
of the monohydrides OH and AgH"^. The OH molecule 
is chosen because precise experimental data [5, 6, 7] as 
well as theoretical data [8] are available. To our knowl- 
edge no experimental data exist for the AgH"^ molecule. 
However due to the fact that the hydrogen atom is 
slightly polarized by the Ag"^ core, the bond is weak and 
an accurate determination of binding energy and bond 
length requires high precision of the theoretical method. 
Our AgH^ result is compared to the pseudo-potential 
calculation of Ref. [10]. 



2. Computational Details and Results 

The first step of the calculation is related to comput- 
ing the atomic wave functions. The HF functions are 
obtained with the one-electron center of gravity approx- 
imation and calculated for each atom of the molecule. 
Instead of taking large-radius unoccupied orbitals in 
building CI wave functions, we have added the Sturm's 
orbitals. They are localized and create a nearly complete 
basis set. A multi-configuration CI technique is em- 
ployed to calculate correlation corrections. For the OH 
molecule two CI basis sets were considered. The first 
has 38 configurations. The configurations are built 
from the Ip"^ valence shell of oxygen, with additional 
single and double excitations of the 2p electrons into the 
3^, 3/7, 3d, 4^, 4/7, and 4d orbitals, and the 1^ orbital of 
hydrogen with excitation in the 2s, 2/7 , 3^ , 3/7 , and 3d 
orbitals. The excited states in both atoms are described 
by Sturm's wave functions. The ls^2s^ shells of the O 
atom are treated as a core. The second basis has 262 
configurations. The CI configurations include single- 
excitations of 1^^, 2s^, and 2/7 "^ shells as well as double- 
excitation of the 2p^ shell where again these excitations 
are represented in terms of Sturm's orbitals. Moreover 
both basis sets have configurations that represent the 
ionic states 0"^H~, and 0~H^. The inclusion of excitation 
from the 1^ and 2s orbitals of the oxygen atom in the 
second basis set leads to a better treatment of the 
core polarization. Table 1 shows that this substantially 
improves the agreement with experiment. The table 
includes the spectroscopical constants of the ground 
state of OH. Our best values disagree with experimental 
values by 0.2 % for the equilibrium internuclear dis- 
tance 7?e, 2.7 % for the vibrational constant We, 0.4 % for 
the rotational constant ^e, and 4 % for the dissociation 
energy D^. We want to point out that even with a very 
small basis set (38 configurations) we reached good 
agreement with experiment for the R^ and ^e para- 
meters. Another theoretical result [8] obtained by a 



Table 1. Constants of the X^II OH molecule (1 A = 0.1 nm; the energy equivalent 1 cm"^ is 29.979 245 8 GHz; HF is the Hartree-Fock method) 





Basis set 












(number of conf.) 


R. (A) 


We (cm ^) 


B (cm-^) 


/)c(eV) 


This work 


38 


0.972 


7297.25 


18.782 


-3.1701 




262 


0.967 


3638.39 


18.834 


^.4302 


HF+^-functional [8] 




0.947 






^.4082 


Experiment [7] 




0.969 


3737.76 


18.910 


^.6259 
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Hartree-Fock method using a A' -functional treatment of 
correlation effects has a larger difference from the 
experimental values than ours. The origin of the dis- 
agreement in the dissociation energy of the VB and HF 
method seems to be related to the fundamental differ- 
ence in constructing the molecular wave function. In the 
VB method the molecular wave function is constructed 
directly from many-electron atomic determinants with 
singly occupied atomic orbitals. For large internuclear 
separation this molecular wave function has the correct 
asymptotic behaviour and reduces to a product of the 
wave functions of the individual atoms. The HF method 
makes use of molecular orbital (MO) wave functions, 
which contain two electrons at all internuclear separa- 
tions. In general, the one-determinant form of the HF 
method is a very poor approximation for large inter- 
nuclear distances. As it is described in Ref. [9] the 
molecular wave function in Ref. [8] is introduced in a 
two-determinant form, which enables them to improve 
asymptotic behavior of the energy. However, this 
modification is still insufficient to get a fair agreement 
of Z)e with the experimental value [7]. 

The molecular constant, /?e, from Ref. [8] and our 
value (Table 1), lie on the opposite sides of the experi- 
mental value [7]. This can be related to the difference in 
the treatment of correlation in the two approaches. In 
our study we apply Sturm's wave functions in the CI 
expansion. This is a very efficient treatment of correla- 
tion effects in atoms and molecules, because the purely 
discrete spectrum of Sturm's eigenvalues contains the 
continuum spectrum as well and hence the CI expansion 
converges quickly. On the other hand, the ^-functional 
technique of Ref. [8], makes use of an effective potential 
obtained by modifying the exchange operator in the 
Hamiltonian. This modification does not effect the 
closed shells of the interacting atoms therefore the 
polarization and correlation effects of the core are not 
taken into account. Finaly, we conclude that the use of 
the Sturm's CI expansions for all orbitals of the OH 
molecule gives us an advantage in calculating the 
equilibrium internuclear distance R^. 



The CI treatment of the AgH^ molecule excludes all 
closed orbitals of the Ag^ ion with the exception of the 
Ap^Ad^^ shells. Single- and double-excitations of the Ap^ 
and Ad^^ shells are considered. Additional Sturm's 
functions for 5^, 5/?, 5d, 6s, and 6p orbitals were con- 
structed for the CI expansion. The number of configura- 
tions is 150. Our calculated values for the molecular 
constants of AgH^ are given in Table 2. The other 
theoretical results [10] are obtained by using a DF pseu- 
dopotential and a density functional approach to de- 
scribe the valence correlation and polarization energies 
of the molecule. Due to the fact that AgH^ is a "single" 
electron molecule, the correlation effects are not very 
important. However, the presence of the hydrogen atom 
polarizes the Ag^ ion orbitals and the CI treatment is 
still required. The comparison with other theoretical 
results [10] shows that we agree very well in the calcu- 
lation of dissociation energies and mostly disagree in the 
computation of the rotational constant We- 



3. Conclusions 

In this paper we have shown that a quite accurate 
determination of molecular constants of molecules with 
open shells using our version of the valence-bond 
method is possible. We treated O and Ag as many-elec- 
tron atoms, where core polarization must be properly 
accounted for. The disagreement with experimental 
values varies from 0.2 % to 4 % for the molecular 
parameters of OH. In addition the calculation of OH 
with a small number of configurations (38), which 
requires a few seconds of computer time for each inter- 
nuclear separation, has less than a 1 % discrepancy for 
the /?e and B^ parameters. The difference between our 
values and other theoretical values for the AgH"^ 
molecule does not exceed 7 %, except for We, which also 
serves as a further test of our method. 
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Table 2. Constants of the X^X AgA"^ molecule (DF is the Dirac-Fock method) 



Basis set 
(number of conf.) 



R. (A) 



We (cm ) 



B (cm-^ 



D. (eV) 



This work 
DF [10] 
DF+core [10] 



150 



2.19 


366.14 


2.37 


515 


2.41 


519 



3.5172 



-0-145 

-0-15 

-0-14 
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